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The higher dimensional Weyl curvature induces on the brane a new source of gravity. This Weyl 
fluid of geometrical origin (reducing in the spherically symmetric, static configuration to a dark 
radiation and dark pressure) modifies space-time geometry around galaxies and has been shown to 
explain the flatness of galactic rotation curves. Independent observations for discerning between 
the Weyl fluid and other dark matter models are necessary. Gravitational lensing could provide 
such a test. Therefore we study null geodesies and weak gravitational lensing in the dark radiation 
dominated region of galaxies in a class of spherically symmetric brane- world metrics. We find 
that the lensing profile in the brane-world scenario is distinguishable from dark matter lensing, 
despite both the brane-world scenario and dark matter models fitting the rotation curve data. In 
particular, in the asymptotic regions light deflection is 18% enhanced as compared to dark matter 
halo predictions. For a linear equation of state of the Weyl fluid we further find a critical radius, 
below which brane-world effects reduce, while above it they amplify light deflection. This is in 
contrast to any dark matter model, the addition of which always increases the deflection angle. 
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I. INTRODUCTION 



The idea of embedding our Universe in a higher dimensional non-compactified space-time has attracted considerable 
interest in the last decade, due to the proposal by Randall and Sundrum [l| that our four-dimensional (4d) space- 
" time could be a three-brane embedded in a 5-dimensional (5d) space-time (the bulk). According to the brane-world 
\Q , scenario, the physical fields (electromagnetic, Yang-Mills etc.) observed in our 4d Universe are confined to the three- 
t-H 1 brane. Only gravity can freely propagate both on the brane and in the bulk, with the gravitational self-couplings not 
C*~) \ significantly modified. The model allows for a large, or even infinite non-compact extra dimension, in the simplest 
£ — , ■ case the brane being identified to a domain wall in a 5d anti-de Sitter space-time. Even with the fifth dimension 
uncompactified, standard 4d gravity can be reproduced on the brane in a certain limit. For a review of the dynamics 
CNI ■ and geometry of brane universes, see e.g. 0. 

At very high energies, in the presence of large 5d curvatures, significant deviations from the standard Einstein 
theory could occur in brane-world models, due to the nonstandard model 5d fields, possible asymmetric embeddings 
of the brane into the bulk or a variable brane tension [H, H|. At the electro- weak scale of about 1 TeV gravity is 
largely modified. The cosmological and astrophysical implications of the brane-world theories have been extensively 
investigated in the physical literature 043 • 

On the vacuum brane the gravitational field equations depend on the generally unknown brane projections of 
the Weyl curvature of the bulk, generating non-local brane stresses, which in a spherically symmetric setup can be 
expressed in terms of two functions, the dark radiation U and the dark pressure P, respectively 0, |H,Q. On a vacuum 
brane several classes of spherically symmetric solutions of the static gravitational field equations have been found in 

Dark matter is one of the central issues in modern astrophysics (see [HI for an extensive review of the recent results 
of the search for dark matter). The necessity of considering the existence of dark matter at a galactic and extra- 
galactic scale is required by two fundamental observational evidences, the behavior of the galactic rotation curves, 
and the mass discrepancy in clusters of galaxies, respectively. On the galactic/intergalactic scale the rotation curves 
of spiral galaxies [14j4l6| provide one of the best evidences showing the problems Newtonian gravity and/or standard 
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general relativity have to face. The rotational velocities of hydrogen clouds in stable circular orbits increase near the 
center of the galaxy, in agreement with the standard gravitational theory, but then remain nearly constant at a value 
of v t goo ~ 200-^ 300 km/s Hence we obtain a mass profile of the form M(r) = rvf goo /G. This result implies that 
the mass within a distance r from the center of the galaxy increases linearly with r, even at large distances where 
very little luminous matter does exist. 

The second important astrophysical evidence for dark matter comes from the study of the clusters of galaxies. 
Generally it is found that the virial mass Myr is considerably greater than the observed baryonic mass Mb, Mvt > 
M B , typical values of M VT /M B being about 20 - 30 [3. 

This behavior of the galactic rotation curves and of the virial mass of galaxy clusters is usually explained by 
postulating the existence of some dark (invisible) matter, distributed in a spherical halo around the galaxies. The 
dark matter is assumed to be a cold, pressure-less medium. Many possible candidates for dark matter have been 
proposed, the most popular ones being the WIMPs (Weakly Interacting Massive Particles) (for a review of the 
particle physics aspects of the dark matter see [l7j])- While extremely small, their interaction cross sections with 
normal baryonic matter, are expected to be non-zero, so that their direct experimental detection may be possible. 

However, up to now no direct evidence or annihilation radiation from dark matter has been observed, and no 
non-gravitational evidence for dark matter does exist. Therefore, it seems that the possibility that Einstein's (and 
the Newtonian) theory of gravity breaks down at the scale of galaxies cannot be excluded a priori. 

In brane- world models the rotational galactic curves can be naturally explained without introducing any supplemen- 
tary hypothesis [ill, E3 • The non-zero contribution of the Weyl tensor from the bulk generates a modified, spherically 
symmetric geometry, in which galaxies are embedded. The dark radiation U and the dark pressure P act similarly 
to a "matter" distribution outside the galaxy. The particles moving in this geometry feel the gravitational effects of 
U, which can be expressed in terms of an equivalent mass (the dark mass) Mjj. The dark mass is linearly increasing 
with the distance, and proportional to the baryonic mass of the galaxy, Mjj{r) ss Ms(r/ro) fill ]. 

Therefore it would be of uttermost importance to have independent tests, which could discern between various dark 
matter models and modified gravity models, which include brane-worlds. 

Gravitational lensing and the study of the light deflection by black holes and galaxies is an important physical effect 
that could provide specific signatures for testing the brane-world models (for a review of the gravitational lensing by 
brane- world black holes see [18|). Observables related to the relativistic images of strong field gravitational lensing 
could in principle be used to distinguish between different brane-world black hole metrics in future observations. 

It is the purpose of the present paper to consider the lensing in the dark radiation dominated region of the brane- 
world models. Physically, this region corresponds to particles gravitating in circular orbits and at constant speed 
around the galactic center [l~4| . 

The galactic rotation curves show a large variety of behaviors, and, in particular, they also depend on the considered 
galaxy type. By analyzing an extended set of spiral galaxy rotation curves it has been proposed that the rotation 
curves of these galaxies can be described by a Universal Rotation Curve (URC) N-body simulations provide 

a universal mass profile in the A Cold Dark Matter (ACDM) cosmological scenario, and, consequently, a universal 
equilibrium circular velocity of the virialized objects, as galaxies. By combining kinematical data of inner galactic 
regions with global observational properties, the URC of disc galaxies and the corresponding mass distribution out to 
their virial radius was obtained. The existence of a Universal Rotation Curve is also consistent with the predictions 
of the brane world models, which imply the existence of of a Universal Weyl Fluid, acting in a uniform way at all 
galactic scales. In the present paper, we will analyze the lensing properties in both the constant velocity region of 
the rotation curves, as well as in the declining or increasing regions, where the behavior of the rotation curves can be 
modeled by a simple power law (l3| . 

By fixing two radius values at the observer and source, respectively, we derive an exact lens equation relating two 
angular variables. This equation allows to obtain all the observationally relevant quantities, like image position, the 
apparent brightness and the image distortions. 

The present paper is organized as follows. We review the basic properties of the brane world models and of the 
Weyl fluid in Section [H] The metric properties valid in the region of constant or power law type tangential velocities 
are also discussed. The mathematical problem of the embedding is considered in Section IIIII Particular solutions of 
the field equations describing the gcocmtric properties in the constant and power law velocity regions are presented in 
Section IIVI The deflection of light and the lensing properties in the brane world models in the Weyl fluid dominated 
regions are obtained in Section [V] In Section I VII we compare the theoretical predictions of the model with the 
observational data. We discuss and conclude our results in Section IVIII 
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II. THE WEYL FLUID IN A SPHERICALLY SYMMETRIC, STATIC VACUUM BRANE 

A. A succinct introduction to brane-worlds 

In the brane-world model standard model fields are confined to the 4d brane, a hypersurface (M, g^) embedded 
(for simplicity, in a Z 2 symmetric way) in a 5d bulk space-time M qab)- We denote the coordinates in the 
bulk and on the brane as X A , A = 0,1, 4 and x M , /i — 0, 1, 2, 3, respectively. The system is characterized by the 
combined action Q 
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where k\ = 87rG 5 is the 5d gravitational constant, is the 5d scalar curvature, A 5 is the 5d cosmological constant. In 
the boundary terms denote the traces of the extrinsic curvatures taken on the two sides of the brane, Lb rane (g a /3, ip) 
is the 4d matter Lagrangian, a generic functional of the brane metric g a p and standard model fields tp, finally A/, is 
the brane tension (chosen here for simplicity a constant). 
The Einstein field equations in the bulk emerge as Q 



^G u = ~A 5 ^ 9IJ + 6(0) 



-A 6 (5) 5 /j + T Ll 
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where gjj is the induced brane metric expressed in bulk coordinates, the Dirac delta-function 6 (0) appears in the 
source terms localized on the brane, including the brane tension and T/j, the energy-momentum tensor of the standard 
model fields, defined in brane coordinates as 
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The brane energy momentum tensor is related to the jump in the extrinsic curvature = Cng^v by the Lancyos 
equation (the second Israel junction condition) 



K liv (x»,y = +0) - K^(x^y = -0) = -k z 5 T f 



{lis 



1 



T 



(0) 



The 4d equations emerge as projections to the brane with the induced metric gu and contractions with the brane 
unit normal n 1 (the brane has the normal form njdx 1 = dy and is chosen at y = 0). The tensorial projection is the 
effective Einstein equation Q 



with the quadratic energy-momentum correction (to be ignored at infrared scales) 

Sp* = 12^^" — 4^» Q ^a + ~^9iiv (3T Q/3 r Q(3 — T 2 ) , 



(7) 



(8) 



and Efj, v a non-local effect from the bulk gravitational field, the electric projection Ejj = Gjajbti n of the bulk Weyl 
tensor Cjajb (with E]j —> E^djS'j as y 0). We denote k\ = 8irG, with G the 4d gravitational constant. The 
4d cosmological constant A, and the 4d coupling constant kn are related to the brane tension as A = fc|(A5 + fc|A 2 /6)/2 
and k\ — fc|Ab/6, respectively. 
The Codazzi equation 
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implies the conservation of the energy-momentum tensor of the matter on the brane, D V T^ U = 0, where D u denotes 
the brane covariant derivative. Moreover, from the contracted Bianchi identities on the brane it follows that the 
projected Weyl tensor obeys the constraint 

D„V = k\D v S^. (10) 
D V E„ V = 0. (11) 



For low density region, we have, 



B. The Weyl fluid 

The 5D vacuum metric containing a spherical symmetric, static, vacuum brane could be written in terms of normal 
coordinates as 

( 5 W = 9AB dx A dx B = -M(r, yfdt 2 + N{r, yfdr 2 + Q(r, y) 2 (dd 2 + sin 2 d(f> 2 ) + dy 2 . (12) 
The corresponding electric part of the Weyl curvature becomes 

r ,( 1 2Q,y,y 2Ny y A 5 Nyy A5 Q y y A5 QjJJJ \ 

+ + ________ (13) 




where we followed the decomposition [22 

N y , y = 6(y)N y + N^ v , (14) 

Q,y,y = fi{y)N,y + Q,y,y, (15) 

into a distributional part and an analytic function (regular) part. The regular part is subject to the constraint (|lll) . 

characterizing the bulk effects on the brane. On the other hand N tVtV and Q, y , y depend on the embedding of the brane, 
they corresponding to the jump in N. y and Q, y , respectively, which arise from the jump in the extrinsic curvature. 
The symmetry properties of also imply the irreducible decomposition with respect to a chosen 4- velocity held 
on the brane [2| 



I ' [ u^u v + ) + P^v + 2V( J1 u u ) 



(16) 



where k = k^/ki, /i M „ — + u^Uu projects orthogonal to u^, the dark radiation U — —k i E^ v u^u v is a scalar, 



Vn = k^h^Eap is a spatial vector and 



p 



h {ll a h v) P - I h^K* 



EaP, (17) 



is a spatial, symmetric and trace-free tensor. 

In a static spherical spacetime V M = 0, and Eq. ([TT]) takes the explicit form Q: 



-D I1 U+-UA I1 + D V P I1U +A V P I1V =Q, (18) 



1 4 
3^+3 

with — u v D v u^ the 4-acceleration. By further assuming spherical symmetry, we may chose A^ — A(r)r^ and 



P^ = P(r) (r»r v - 1 v) . ( 19 ) 

where P(r) is the dark pressure (although the term dark anisotropic stress might be more appropriate) and r M is a 
unit radial vector [8(. 

Characterized by just two functions U and P , the electric part E^ u of the Weyl curvature can be regarded as a 
perfect fluid of geometric origin, the source of gravity on the brane. We call this "fluid" the Weyl fluid. 

The Weyl fluid is a macroscopic fluid with an effective energy density and pressure, similar to any other material 
fluid. However, the energy and the pressure of the fluid depend and are determined by the geometry of the extra 
dimension. Moreover, it satisfies a more complicated continuity equation than usual material fluids, Eq. (|10p . and 
it is traceless. If we observe a Weyl -type fluid behavior in the Universe, the most natural explanation will be the 
brane world effects. Using the hypothesis of the Weyl fluid as a solution for the missing mass problem is a different 
approach than assuming the existence of a new particle, generating a standard material fluid in general relativistic 
framework. 
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C. Dark matter as a bulk effect 



In this subsection we discuss a brane metric properly reproducing the galactic rotation curves. We assume that 
galactic dynamics is governed by the static and spherically symmetric 4d line element 

ds 2 = -e v ^dt 2 + e x{r) dr 2 + rW, (20) 

where dVl 2 — d6 2 + sin 2 9d(j) 2 . Here and </> are spherical coordinates, t € R and r ranges over an open interval 
(/min, '"max) so that < r m j n < r max < co. We also assume that the functions v(r) and A(r) are strictly positive and 
(at least piecewise) differentiable on the interval (r mm , r max ). 

In the metric (|20p the independent components of the effective Einstein equation are 

^ xfl X '^'-± = -Ei, (21) 



r 

/,.' 1 

-A 



^ ( v - + i ) - i = -k, (22) 



e x ( ., v' 2 v 1 - a' 



(23) 

where a prime denotes the derivative with respect to r. Knowing the brane metric, the projected Weyl tensor is 
determined by this system. 

The tangential velocity v t„ of a test particle is measured in terms of proper time, that is, by an observer located at 
the given point, as [20I [2lj 

v 2 g = e-"r 2 n 2 . (24) 



By using the constants of motion, the tangential velocity of a test particle in a stable circular orbit on the brane 
emerges as [2(J [H| : 



2 



rv 



(25) 



"tg- 2 ■ 

In a power law tangential velocity limiting case fl3j |. 

v tg (r)^v c r^ (26) 

where v c and C are constants. 

This power law velocity profile could be the limiting case of various velocity profile. For example, such a behavior 
is suggested by the Universal Rotation Curve model [lM [l|| . The corresponding metric coefficient is given by 

^ + „ 0) for C + , 

> (27) 
In (-^-J , for C = 0., 

where i/ and R^, are arbitrary constants of integration. 

In particular, constant tangential velocity limit immediately allows us to find the metric tensor component e v in 
the flat rotation curves region on the brane as 



r 



R 

Since is traceless, summing all Eq. ([2"Tj) -Eq. (|2"5)l could determine 



(28) 



e" A «l-<. (29) 



Eqs. (l28l) .and (l29l) give the galactic metric in the constant velocity region on the brane 
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The brane metric obtained phenomenologically by the requirement that the required tangential velocities are realized 
in the model could furter be specified by assuming a particular embedding of the brane in the bulk, equivalent with 
the equation of state of the Weyl fluid 

P= {a -2)U-^. (30) 

Here 3a b = 64tt 2 G 2 and B is an arbitrary constant of integration. It has been shown, that this choice of the equation 
of state leads to a tangential velocity compatible with the observations on rotation curves [TT1 . [23| . An explicit proof 
based on the 3+1+1 decomposed covariant bulk dynamics [3| has shown that the corresponding five dimensional 
embedding does exists |23| . 

In the next Section we will discuss the embedding of this brane metric into 5d. 



III. THE EMBEDDING PROBLEM 



In the previous Section we have used the observed characteristics of the rotation curves to fix e v , and the traceless 
property of the electric part of the Weyl tensor to fix e A . We need to prove that our choices are compatible with an 
integral bulk solution, i.e. the derived phenomenological metric can be realized on a brane embedded in some bulk 
geometry satisfying the 5d Einstein equations ((4]). 



A. The Campbell-Magaard theorem 



The Campbell-Magaard theorem states that it is possible to embed a manifold with arbitrary geometry into a 
Ricci-flat manifold with one extra dimension [24], [25[ . The use of the Campbell-Magaard theorem for the embedding 
problem in higher dimensional theories has been studied recently (27l - [30l |. In [26[ it was shown that it is possible 
to generalize the Campbell-Magaard theorem to the Randall-Sundrum type brane-world model, as any solution of 
4d projected field equation can be realized as a thin 3-brane in the brane-world. However, if one considers the 5d 
geodesic corrections to particle dynamics, the matter on brane will leak off if the total energy momentum on the 
brane (including the brane tension) does not satisfy the strong energy condition. Although we do not consider higher 
dimensional particle motion, we include a discussion on the embedding based on the Campbell-Magaard theorem. We 
find that a metric reproducing rotation curve dynamics could be realized in the Randall-Sundrum type brane-world, 
i.e. there is a 5d spacetime (M ,gMN) satisfying the field equations such that the brane is a representant of a foliation 
of M with codimcnsion 1, giving the rotation curve dynamics. In the following subsection, we will give the details 
based on the procedure described in [26]: 

1. Write down a brane metric g^ v , extrinsic curvature K^ v and projected Weyl tensor that is consistent with 
the equations for §mn evaluated on the brane. 

2. Show the existence of a unique solution for qmn on the half plane y > with boundary condition given in 1. 

3. Show the existence of another unique solution for equations of guN on half plane y < 0. 

4. The union of the two solution will be our bulk spacetime. 

However, the bulk geometry will not be AdS5, and thus gravity proves to be localized on the brane, with a Newtonian 
limit. This is acceptable, as there is no independent experimental test for gravity and particle physics on the galactic 
scale. In the following subsection we discuss the details of this procedure. 



B. Existence of the bulk solution 



We will make use of Cauchy-Kowalewski theorem for non-linear partial differential equation (PDE) [31j . The 
theorem states that a unique analytic solution exist in the neighborhood of y = for an analytic Cauchy problem. 
The analytic Cauchy problem is defined by a system of PDEs with a Cauchy data. The PDEs have the form 

d™u = F i {x,y,u{x>y)M a) ) (31) 
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where x can be an array, is the derivative of the dependent variable u to order a < m, and F is analytic. The 
Cauchy data are a set of boundary conditions at y = 0, 

d k u = f k ; < k < 77i (32) 

where f k is analytic for every k < m. 

C. The Cauchy problem in the braneworld model 

We are going to translate the braneworld embedding problem into a Cauchy type problem. In the brane world 
model we are considering a weak solution of the Einstein equations, i.e., gMN^x^, y) satisfies the 5D Einstein equation 
for a region y ^ 0, and the brane is a discontinuity on d y gMN(x^ , y) at y = 0. Therefore, we first look at the solution 
of the 5D Einstein equation inj/e [0, +oo). 

On the other hand, using the function-distribution decomposition, Eq. ()14j) . the 5D Einstein equations become 
PDEs without delta functions that are valid on each half plane. 

(5) 2Q,r-N^ r ^ 2QyNy ^ 2Q TT ^ Nyy 



where 
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Q 2 N 2 Q 2 Q 2 { ' 



This system forms the analytic PDEs of the Cauchy problem. 
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D. Analyticity and consistency of Cauchy data 

We need to show that the way we fixed the 4d metric in our previous Section gives consistent Cauchy data. We 
determine e" by the rotation curve and determine e A by TrG M „ — 0. On the other hand the Cauchy data relevant to 
Eqs. ([55 ]l -([5o j) consist of g(x»,0) and 9^(^,0). 

The dynamics inside galaxies, says for r < r c , is not described by our brane metric. The metric specified by 
Eqs. (I2TJ1) . (|2"4"|) and the traceless property of E^ v only applies to the outer region of galaxies, so that they have a 
domain of r £ [r c ,oo). Therefore the physical situation under consideration isolates singularities that may arise in 
the brane metric, i.e., fixing the 4d metric determines g(x fJ ',0), which is analytic on r 6 (r c ,oo). 

On the other hand, the junction condition determines d y g(x fl , +0), which is analytic. However, g(x ,1 ,0) also tells 
about dy y g(x^ 1 , 0) by equating Eq. (fl3]) with Eq. ([21"]) - Eq. (f23]) . But these equations must be consistent with 
Eqs. (|33| -(|36 |) as Eqs. (fT3|) are derived from their geometric projection. One could also verify these by inserting them 
into Eqs. (Pl) - (j5rJT) . 



E. Zi symmetry and AdSs bulk 

We know that a unique solution on the half plane y € [0,+oo) does exists. By putting g(x^ l ,z) = g(x^,—y) into 
Eqs. ([33l) - (j36]) . we obtain the PDEs on the another half plane. g(x^ 7 0) and d y g(x ,i , z = +0) give the Cauchy data for 
the another half plane. The uniqueness of solution, the symmetry of Eqs. and the Cauchy data shows that 

the solution on the another half plane is identical, i.e. g(x^, —y) = g(x ti ,y). 

However, there is no guarantee of an asymptotic AdS5 solution. Such a requirement would over-determine the 
Cauchy problem. This can be seen, for example, from the numerical analysis of the spherically symmetric solution 
[32I ]. A numerical solution for 5d metric that is ADS asymptotically can be obtained without any phenomenological 
assumption on the brane metric. As we have discussed in the Introduction section, experimental confirmation of 
Newtonian gravity on galactic scale is still lacking. 



IV. PARTICULAR SOLUTIONS OF THE FIELD EQUATIONS IN THE WEYL FLUID DOMINATED 

REGION 

In the present Section we present particular solutions of the gravitational field equations on the brane in the 
dominated regions (where we include baryonic matter only as test particles). 

The metric that describes the spacetime in the spherically symmetric galactic halo region, with the baryonic matter 
concentrated in the central part of the galaxy, is the Schwarzschild solution of general relativity. A brane-world 
solution could emerge as a correction to the Schwarzschild metric, with 



e "M = i_l£ +jBt ( r ) (38) 
r 



e A(r) = ^ rr ■ (39) 

Here rs is the Schwarzschild radius. The functions B t {r) and B r (r) should be compatible with a traceless Einstein 
tensor G^ v in order to be realizable in Randall Sundrum brane-world model. The traceless condition gives rise to the 
following equation 

TrG = B' t '(r-r s + rB r ) _ rB' t 2 (rB r + r - r s ) rB' r B' t (r - r s + rB t ) 
r-rs + rBt 2(r-r s + rB r ) 2(r -r s +rB t ) 

The question comes whether it is possible to have a solution for B t , B r consistent with the mass required from 
galactic rotation curves. For this, B t {r) should be related to the tangential velocity profile Vt g (r) of the galaxies 
through 

vl(r) = ^. (41) 

Substituting this equation into the metric (f3"5)) , we obtain the differential equation for B t in terms of the tangential 
velocity profile 

B' t - 2 ^(l- r ^ + B t ) + r 4 = 0. (42) 
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Eliminating B' t from Eq. (|40|) leads to a differential equation for B r , 

(2 + OK + (J + ^ + ^ + ^) B r + 4 a (I - ^) + < - 3) + ^ (l -2) = 0. (43) 

The solution of this system of ordinary differential equation gives the corrected Schwarzschild type solution of the 
Einstein equations. 

A. The metric in the flat rotation curves region 

In the limit of constant tangential velocity, i.e. v tg — > Uqo, Eq. (j42|) can be integrated and leads to the asymptotic 
behavior 

B t (r)^-l + ^ + C t r 2 ^. (44) 
r 

Eq. (|4"3")l can be integrated to give 



B r (r) r [C r - r "~+ 2 



^4 + 1 



rs 



(45) 



with Ct, C r being integration constants. 

The virial mass is derived from the observed acceleration of massive baryonic particles, for example, in X-ray cluster 
images. It will be useful to see the acceleration of massive particles in the weak field, low velocity limit of our metric. 
Massive particles move on brane geodesies 

fft T P rl r P A T <J 

— — + r M ^L^L - o (46) 

dr 2 pa dr dr ' { ' 

with t the proper time. At small velocities dx^/dr rj (1, 0, 0, 0) and t s» r hold [33j, and thus 

d 2 x» 



dt 2 



K , (47) 



a=^(l- r i + B r )(^ + rBl) . (48) 

Inserting the expressions for B t , B r from Eqs. (|44l) . (|45[) and keeping the leading order terms in the small parameters 
rs/r and gives 

v 2 

a « -2^ (49) 
r 

The acceleration derived from this brane-world metric agrees with the Newtonian gravitational acceleration of dark 
matter halo with inverse square law radial dependence of the gravitational force. 

B. The unified Schwarzschild - constant velocity metric 

We consider here a tangential velocity arising from a baryonic central mass and some unknown mass as 

«*(r) a = y+i&- (50) 
This tangential velocity profile could be described by the metric 

2uf 



e 



^ = (~) ^fl-?), (51) 



r, 



eA(r) ~ r^HW' (52) 
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where r c is a characteristic radius of an individual galaxy. Using the traceless property of dark radiation we obtain 
the differential equation for B r , 



vlr + 2r- Irs ) K + 2 ( 1 + *4 - — ) B r + 2t& = 0, 



(53) 



where we have used relevant approximations valid for the galactic motion. B r can be obtained by power series. 
We also investigate the acceleration of slowly moving massive particle by estimating the Christoffel symbol, 



r - i [ r - 



2oZ 



r 



2v °c ( 1 j + — 



1 / r 
~2r \7r. 



2ri 



^- + 7)«^ + r- (54) 

r / r 2r 



The acceleration of the particles agree with the Newtonian gravitation acceleration produced by a central mass and 
a dark matter halo. 



C. The metric in the power law velocity profile region 



From Eq. (f!?T|) with £ ^ 0, we know directly one metric coefficient 



In the case Q < 0, we can take the large radius approximation, 



? « A 1 - 



with 77 = This is equivalent to a Schwarzschild correction in the asymptotic region, 

Bt(r)->-l + — ''" 



2), 



(55) 



(56) 



(57) 



The other correction _B t can be obtained in large radius limit from Eq. (|43p , which can be simplified in this approxi- 
mation as 



This metric coefficient gives the radial acceleration as 

IT, « Av1r-^ +1 \ 



(58) 



(59) 



If A = 1, the acceleration is consistent with the acceleration produced by the dark matter density profile suggested 
by the power law velocity profile 



(60) 



D. Solutions with dark pressure equation of state 

For the equation of state (|30]) of the Weyl fluid, the solution of the effective Einstein equation can be written as 

e"*M = 1 - ^ - GjV% (r) , (61) 



exp 



2a - 3 

7 

a — 3 



(62) 
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where 



GM v (r) « v r + d + C 2 r^ - 2GM , (63) 



and Cb,Cu,Ci, C2 and M are also arbitrary constants of integration. The approximate equalities indicate that in this 
approximation baryonic matter (which is subdominant), was dropped. Finally the constants vq and the asymptotic 
tangential velocity v tgoc are given in terms of the previously introduced constants as 

B(B-Z) . . 

U0 = a(2B-3) + 9 ' (64) 



v t goo = 3 ( av o ~ B) . (65) 

This solution has been presented in detail in 23] . 

In the weak field regime we have further simplifications, as has been discussed in (23[. A careful post-Newtonian 
counting first gives B <C a — 2, then 

V ~ Vtgoo ~ ~^ . ( 66 ) 

with either a < 3/2 and B < or a > 3 and i? > 0. 

The constant C„ in Eq. (|6lT) can be immediately absorbed in the time coordinate. We introduce the following 
notations 



B 



2a — 3 ' ^ a — 3 ' 
C = Car?" 1 , GM Q = ^ , (67) 

and two mass type constant Mq and M& given by 

2GMi = C* 6 + 2G (M - M) , (68) 
The metric functions with the speed of light c reintroduced can be approximated as 

e-W«l- 7 -^-c(iy~\ (69) 



e" v '» I — I exp 



a f(r) 



2GM 2C 



1 



(70) 



In general Afi and Mq are two independent parameters, however Mi does not affect the rotation curves and should 
be determined by other observation, like gravitational lensing. In the rotation curve studies Mi = Mq has been 
chosen. The rotation curves explained in terms of the Weyl fluid were characterized by three universal dimensionless 
constants a, 7, C and a mass-type constant Mq. The constant r c was introduced in order to have C dimensionless. 
In LSB galaxies the visible matter can be considered to be concentrated inside a constant density core, as shown in 



34] . In a stellar and gas free model, the matching of the outside brane- world solution to the inside Newtonian region 
suggests that C = —7 and below r c visible matter dominates. 

The constants a, 7 obey either a<0or0<a<l (the value a = is excluded, as it would correspond to the 
unphysical values a — > ±00. In both cases 7 is a small positive number < 7 <C 1 (with the exception of a ~ 3, 
translating to a ~ 1, when theoretically 7 can be an arbitrary positive number). The choice from the analysis of the 
Low Surface Brightness (LSB) galaxy rotation curves [23| and the fitting from LSB galaxies suggest that 7 < 1. The 
values of the parameters a, 7, r c and Mi are presented in Table HI 
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Galaxy 


Mo 


r c 


a 


7 


Xmin 




M Q 


kpc 








DDO 189 


4.05x10 s 


1.25 


0.3 


6.43x10" 


8 


0.742 


NGC 2366 


1.05xl0 9 


1.47 


0.8 


1.12x10" 


7 


2.538 


NGC 3274 


4.38xl0 8 


0.69 


-0.4 


6.73x10" 


8 


18.099 


NGC 4395 


2.37xl0 8 


0.71 


0.9 


3.43x10" 


7 


27.98 


NGC 4455 


2.26xl0 8 


1.03 


0.9 


2.72x10" 


7 


7.129 


NGC 5023 


2.69xl0 8 


0.74 


0.9 


4.53x10" 


7 


10.614 


UGC 10310 


1.28xl0 9 


2.6 


0.4 


1.12x10" 


7 


0.729 


UGC 1230 


3.87xl0 9 


3.22 


-1.7 


1.12x10" 


7 


0.539 


UGC 3137 


5.32xl0 9 


3.87 


-0.5 


1.23x10" 


7 


4.877 



TABLE I: The best fit parameters (Mo, r c , a, 7) of the sample of 9 LSB galaxies given by [23I ] 

V. GRAVITATIONAL LENSING IN THE WEYL FLUID DOMINATED REGION 

A. Gravitational lensing in brane-world model 

In the present Section we will consider the lensing properties in the dark radiation dominated region of a brane- 
world galaxy, outside the galactic baryonic matter distribution, in the various approximations discussed in Section 
IV. 

We consider the lensing geometry represented in Fig. [TJ with observer O, lens L, source S, image /, and consider 
the deflection of light occurring at a point A. Here 9 is the apparent angle while j3 is the real, undeflected angle of 
the source, and 5 is the deflection angle of the light ray. We denote by D/ s the lens to source distance, and by Di the 
lens distance, measured from the observer. 




FIG. 1: The lensing geometry. 
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The deflection angle 8 is obtained by comparing the two asymptotic directions of the null geodesies [3£ 



S(r B 



Xdr — 7r, 



l(r) = 



9rr (r) 



r V \{-9tt (nnin)) / {-git {r))\ (r/r min ) 2 - 1 ' 



where r m ; n is the distance of minimum approach given by 

<//• 



dep 



(r m in) = 0. 



(71) 



(72) 



(73) 



B. Dark matter lensing in standard general relativity 

In the weak field limit we can use the superposition principle, and therefore the deflection angle for light in a dark 
matter model is given by 



AG 



[M + Mdark(rmin)] , 



(74) 



where M^^r) is the effective mass of dark matter inside radius r. The existence of dark matter always enhances 
the deflection angle. In the dark matter model, we can separate the contribution of the baryonic matter and of the 
dark matter as 



where 



3dark 



5 = 5 S + <S dark , 
4GM dark (r 

min ) 



^min 



For example, if we consider the pseudo-isothermal halo model, 



Piso{r) = po 



r 



(75) 
(76) 

(77) 



where p is the core halo density, and is the radius of core halo, the asymptotic velocity of this model is a flat 
rotation curve, or Eq. (|50p . which is obtained from the unified Schwarzschild - constant velocity metric. 



C. Deflection angle in the flat rotation curves region 



The deflection angle for the brane- world corrections (j45|) and (|44|) can be calculated from Eq. (|71|) . The integration 
constant Ct does not enter the expression of 5, and for the simplest case we also set C r = 0, such that the expansion 
of 5 for small dimensionless parameters r$ /r m - m and is given by 



6 = 2 



1 



^ ^*min 



[i- 



— v 



2(7 



. ) 



dr — 7T. 



(78) 



We see from Eq. (|78p that the first order term in rs/r m ; n vanishes, i.e., the lensing effect of the central masses are 
suppressed by brane- world effects. The integral can be calculated to give for 5 the expression 



o = = constant. 



(79) 



A constant deflection angle can be interpreted as resulting from a galaxy that contains invisible mass with density 
profile 
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The 1/r 2 profile agrees with what standard dark matter models predict. However, if we consider the deflection angle 
as a function of in the dark matter picture, a 1/r 2 density halo gives 

<W = 4<4, (81) 

which is 18% different from the prediction of the brane-world model. Thus brane-world models predict a systematic 
deviation of the dark matter core density deduced from lensing studies as compared to the virial mass obtained from 
the study of the rotation curves. 



D. Deflection angle in the unified Schwarzschild - constant velocity metric 

For the brane geometry described by the metric Eq. (|51j) . the lensing deflection angle can be calculated similarly. 
The post Newtonian expansion of the deflection angle is 




which can be numerically integrated to give 

,5 = 2.000— +4.712*4. (83) 

' min 

The first term is 4GM/c 2 r m ; n , the standard GR result for the deflection angle caused by the central mass. The 
second term is a constant proportional to v 2 ^ , which is equivalent to superposing on the central core an inverse square 
density dark matter halo. However the proportionality constant 4.712, which agrees with the value given in Eq. (|79[) 
is different from that arising in the standard dark matter picture, and as such it could test the model through the 
correlation of the rotation curve data with lensing data. 



E. Deflection angle for the power law velocity profile 

For rotation curves with vt%{r) — v c r~ ri , with rj > 0, the metric coefficient in the braneworld model is given in 
Section IIV CI The term of the order of v 2 in the series expansion of the deflection angle is 

<rj ,22, r a(l-u 2 )-u 2 {l-u 2 «) 

0(r min ) - U c r min _ 1)3/2 ( u + l)3/2 u 2a+l aU ' V* 1 ' 

By taking rj — 0.3 [l3| we obtain 

S(r min ) = 0.345i; 2 r^f n . (85) 

If the density and velocity profile are derived from the gravitational lensing by using general relativity, the velocity 
profile imply a deflection angle 

S(r min ) = Av 2 c r°J n . (86) 

The profile has the same form for both theories, with different characteristic constants. This is a very important 
consistency check of the missing mass model. If dark matter is the correct description of the missing mass problem, 
the velocity profiles and the lensing profiles need to be consistent. 



F. Lensing for models with a Weyl fluid equation of state 



/ \ 27 

In Eqs. (|69l) and (|7D)) . the factor I — ) can be approximated by 1 for 7 -C 1. As Mi in general affects the lensing, 

we assume Mi 7^ Mo, and we set Mi = Mo + Mb, where M& is a brane-world correction to the Schwarzschild central 
mass in g rr . When M& = we reobtain the results for the rotation curves in 23]. Taking into account that the 
exponent in Eq. (|70p is small, i.e., 



2GM 2C 



«1, 



(87) 
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and Eqs. (|69[) and (|70[) can be approximated by 



5 w(r) 



2GM 2 7 



1 — a \ r. 



and 



_ AM 2GAf 2GM fc 



a-l 



(89) 



respectively, where for simplicity in the following we omit the subscript of Mq. Mb cannot be fixed by rotation curves 
alone, and hence it should be constrained by lensing observations. Notice that Eqs. (1551) and (1591 are not valid for 
r — > oo. Since Eqs. (|70|) and (l69|) apply to empty space with the boundary condition given by Eq. (|65|) . the metric 
given by Eqs. (j55J and (1891) should be valid for r c < r < r^, where r m corresponds to the maximum spatial extension 
for which the geometry that describes the rotation curves can be probed. 

Expanding the integrand in Eq. (|T1 [) to first order in both the dimensionless parameter e = GAf/c 2 r m ; n and the 
dimensionless brane-world parameter Eb — GMb/c Vmin, and 7, we find 



X(r) = 2o(r) + sl s (r) + e h 1w\{r) + jlw2 (r) 



where 



Mr) 



Mr) 



r V (r/r min ) 2 - 1 ' 



1 + r/r min + (r/r min ) 5 



1 + r/r 



Mr), 



I\vi(r) = -^-M r )' 
r 



and 



lw2(r) = 



Mr) 



a-l 



(1 - a) [(r/r min ) 2 - 1] 



(90) 

(91) 

(92) 
(93) 

(94) 
(95) 



respectively. From the above series expansion it follows that the deflection angle can be decomposed into the 
Schwarzschild contribution and the brane-world contribution, 



5 = 6s + Sw- 

5s can be integrated by means of the substitution u — r/r min (then also du/u = dr/r), so that 

4GM 



Ss = 



POO 

2/ [l (r)+el s (r)]dr-7r = 

J ^min 



2 : 
C ^*min 



5w can be written as 



(96) 



(97) 



Sw 



2GMb 

[eblwi (r) + ^1wi{r)\ dr = — h 7 



a-l 



where 



9(a) 



du 



u V u 2 — 1 



(1 - a)u a - 1 + (1 + a)u a+1 - 2u 2 



(1 - a) {u 2 - 1) 



(98) 



(99) 



is a function depending on a only. 
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G. Lensing deflection angle 

The contribution of the brane Weyl curvature in Eqs. (|88| and (|89[) (the term proportional to 7) makes (—gtt) to 
increase, while it makes g rr to decrease. Therefore the contribution of the terms containing 7 to I can be of either 
sign. This result can be verified by using the series expansion in Eq. (|98p . The positive or negative sign of 6^ depends 
on the parameters of the system. The effect of dark matter to the lensing is mimicked by the brane- world effects when 
5w is positive. Fig. [2] shows how 8y/ can be positive or negative in the parameter space. 

The existence of a negative Sw in the parameter space is a unique feature of the brane- world models as compared to 
the conventional dark matter models. Standard dark matter models assume an invisible mass distribution surrounding 
the galaxy, with the dark matter particles feeling the gravitational interaction, but interacting weakly with baryonic 
matter in other interactions. The existence of dark matter should always enhance the gravity, and thus increase the 
deflection angle relative to a "visible matter" lens. 



VI. GRAVITATIONAL LENSING AND ROTATION CURVES FOR LSB GALAXIES 

The braneworld Weyl fluid model defined by Eq. (|30[) can be discriminated with respect to standard dark matter 
models by lensing studies. In (23[, the rotation curves of the LSB galaxies were fitted with the Weyl fluid model. We 
would like to compare standard dark matter lensing with braneworld lensing by using the data obtained from the 
corresponding rotation curves fitting in dark matter and braneworld models. 



A. The rotation curves fitting 

The rotation curves of the LSB galaxies used in [23[ have been obtained in |34| . As an example of the standard dark 
matter approach we consider the pseudo-isothermal halo model discussed in [34| . The values of the parameters po and 
ru obtained with the minimum-disk assumption for our sample of 9 LSB galaxies are presented in [34| . However, for 
the comparison of the observational data with the brane world model predictions we only need the qualitative result 
that dark matter always gives positive deflection angle, i.e. it increases the deflection caused by the central mass. 

Fig. [3] represent the lensing profile of the chosen sample of galaxies. The shape of the plots are similar for all 
galaxies in the sample, and we pick po = 10~ 17 g cm -3 and r c = 1 kpc for illustration. 
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10 



oe=0.6 ■ 



a=0.5 ■ 



a=0.1 



a=-0.5 



FIG. 2: The plot shows the evolution of the braneworld contribution Sw to the deflection angle, as function of r m m/r c for 
various values of the braneworld parameter a. The contribution Sw is positive for orbits with the closest approach r m i n above 
the critical radius r^ n (defined by the intersection with the horizontal axis). 
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FIG. 3: The standard dark matter model deflection angle <5dark at different impacts r m i n predicted by the pseudo-isothermal 
halo model, with po = 10 _17 g cm -3 and r c = 1 kpc. We only need the result that the deflection angle is always positive. 



B. Brane-world deflection 



We first define the critical radius of approach by Sw(r m in) = 0. The lensing by Weyl fluid could be sufficiently 
different from that produced by dark matter if lensing occurs inside , as can be seen from comparing Fig. [5] with 
Fig.E 

is given by the equation 



2r b 

r 1 



(100) 



The location of depends on rj, a, r c and 7. From these parameters, only rj cannot be fixed by the rotation 
curve analysis. If is located inside the range of validity of the metrics given by Eq. (|88l) and Eq. there can 
be hope for detecting it. 

In order to investigate the possibility of observing lensing with r m i n < r^*, we pick galaxies DD0189, NGC3274, 
NGC2366 and NGC4455 for analysis. There are two possibilities that could arise from the sample presented in Table 
HI We use the Weyl fluid parameters fixed by rotation curve studies (23[, and plot Sw versus r min for different (in 
units of 10~ 8 r c ). The solid line represents the case n, = 0. There are two cases, illustrated on Figs. 4 and 5. 

In this well hidden inside the galaxy, and therefore gravity is not weakened, but rather 



a. 



r crit <; 



Tr 



y,crit 

' min pa ' c AAi u "" u min 

enhanced anywhere in the asymptotic region. Brane-world contributions and dark matter give similar predictions 

r crit 
min 



b. r^ln > r c In this case could be detected by astronomical observations. The lack of observations of such 



lensing effects may possibly suggests that the choice of the equation of state for the Weyl fluid given by Eq. (|3T))) may 
not be appropriate for the respective galaxies. 



VII. DISCUSSIONS AND FINAL REMARKS 



Dark matter searches have been unsuccessful so far. Therefore attempts to explain galactic rotation curves and 
galactic cluster dynamics by modifying general relativity at large scale cannot be excluded a priori. The spherically 
symmetric brane-world model offers a solution for the missing mass problem without the need to introduce unidentified 
forms of matter. Whether such a model is viable, can be decided by working out a number of its predictions and 
confronting the theoretical predictions with observations. A consistency check on rotation curves and lensing data, 
proposed in the present paper, is therefore relevant for this purpose. 

This paper has shown that specifying an induced metric on the brane yields a unique solution of the bulk in this 
particular problem. This happens because we only consider modification of gravity in the dark radiation dominant 
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FIG. 4: The lensing profile of DDO 189 with different r;, (in unit of 10 8 r c ). Rotation curve data extended to 9kpc, which is 
similar to 7r c . 



region, i.e. the region in which induced metric is analytic. The schemes described in the Campbcll-Magaard and 
Cauchy-Kowalewski theorems allow to locally embed the brane into a region of 5d space-time with negative cosmo- 
logical constant. Therefore we could employ the rotation curve data of galaxies to fix the brane metric. However an 
asymptotic AdS5 bulk geometry is unlikely. 

The corrected metric arising in the brane- world model could explain observations related to the motion of massive 
particles in stable circular orbits around galaxies. Traceless dark radiation could indeed behave like dark matter. The 
model could in principle also explain the missing mass in gravitational lensing without assuming dark matter. 

However, by correlating the rotation curve data and gravitational lensing, the brane- world models show a distinctive 
feature over standard dark matter models. The proportionality constant between the limiting constant tangential 
velocity square and the deflection angle are biased by 18%. This effect could be however blurred by the error margins 
in the rotation curve data. The bias is valid for all galaxies with a flat rotation curve. Although the effect is hard to 
detect in individual galaxies, one can still do a survey on rotation curve-lensing correlation, and investigate the best 
fitting result. 

The unified Schwarzschild - constant velocity metric represents a different way of fixing the brane-world metric 
by using the properties of the rotation curves. It presents another expression for the tangential velocity limit of the 
rotation curves. The present study suggests that one could use simple analytical functions to fit the rotation curves. 
It consolidates the proposal of using observational data on the rotation curves as a probe of the existence of the 
extra dimensions. The predictions of the lensing in the presence of extra-dimensional effects are consistent with the 
limiting constant tangential velocity calculations. Particle motion can be described by introducing an effective mass 
component, but there is some disagreement in the rotation curve and lensing correlation. 

In the model with dark radiation equation of state, we have found that there are certain parameters of the model that 
mimic the observed lensing profiles of the dark matter halos. Galaxies DD0189 and NGC3274 give enhanced gravity in 
all cases. Some numerical values of the parameters in this model however could produce a negative contribution to the 
deflection of light, like, for example, in the case of the galaxies NGC2366 and NGC4455. The negative contribution 
to deflection is due to the repulsive effect of the brane-world gravity; and we consider this an indication that the 
dark radiation equation of state may be inappropriate. The fact that we always found gravitational enhancement in 
galaxies may rule out the choice of the linear equation of state of the Weyl fluid for these galaxies. 

On the other hand, we found that a traceless energy momentum tensor on the brane could explain astronomical 
observations independently of any equation of state of the dark radiation, which could solve the missing mass problem. 
Such an effective fluid could have pressure, but does not have any free streaming. Conventional constraints on dark 
matter that assume that dark matter has to be cold due to the free streaming processes in structure formation history 
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FIG. 5: The lensing profile of NGC 3274 with different r b (in unit of l(T 8 r c ). Rotation curve data extended to 9kpc, which is 
similar to 7r c . 



may not apply to the case of the effective Weyl fluid. 

Obtaining the full extra dimensional metric is among the biggest challenges in brane-world models. Further con- 
straining the equation of state for the dark radiation could also help to better understand both the extra dimensional 
features of brane- worlds, and also the properties of dark matter. 
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FIG. 6: The lensing profile of NGC 2366 with different rs (in unit of 10 8 r c ). Rotation curve data extended to 6kpc, which is 
similar to 4r c . 
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FIG. 7: The lensing profile of NGC 2366 with different r s (in unit of l(T 8 r c ). Rotation curve data extended to 6 kpc, which 
is similar to 4r c . 
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